7个种源蒙古莸叶片解剖结构及地理环境数据集的
内容与研发

纪若璇,于 笑,常 远,沈 超,白雪卡,夏新莉,尹伟伦,刘 超*

北京林业大学生物科学与技术学院林木育种国家工程实验室,北京 100083

摘 要:长期受到生长环境影响而形成的遗传变异对植物生长发育有着显著的影响。叶片是植物对环境变化最敏感的器官,了解叶片解剖结构在不同环境中产生的适应性变异是探索植物对环境适应的基础。同质园试验是研究遗传与环境因素对植物生长代谢等影响的一种有效方法,该研究利用同质园试验排除了环境梯度的影响,通过常规石蜡切片、多重比较、相关性分析、一般线性模型分析等方法,对内蒙古阿巴嘎旗、阿拉善左旗、东乌珠穆沁旗、鄂托克旗蒙西镇、凉城县、河北康保县、陕西神木市7个不同种源地的蒙古莸(Caryopteris mongholica)叶片解剖结构及其影响因素进行了定量比较。同质园试验说明,种源地气候差异驱动的遗传变异是引起不同种源叶片解剖结构差异的重要因素。该数据集包括:(1)种源地地理位置数据;(2)蒙古莸叶片横切面图片;(3)蒙古莸叶片解剖结构与地理环境影响因子分析数据。数据集存储为.shp、.kmz、.xls、.jpg和.doc格式,由14个数据文件组成,数据量为5.84 MB(压缩为1个文件5.75 MB)。

关键词:蒙古莸;同质园;叶片;解剖结构;种源;内蒙古

DOI: https://doi.org/10.3974/geodp.2021.01.13

数据可用性声明:本文关联实体数据集已在《全球变化数据仓储电子杂志(中英文)》出版,可获取:
https://doi.org/10.3974/geodb.2021.01.03.V1.

1 前言

植物在进化过程中,受到外界温度、水分、光照等环境因素的长期作用而影响其建成逐步形成适应环境的表型和遗传特征。即使是同一物种,长期生长在不同的地理环境中,也会产生不同程度的基因分化,以适应各自不同的生态位,从而形成特定的种源特性[1–4]。叶片作为植物暴露在生活环境中面积最大的器官,其外部形态特征、内部解剖结构等最能体现植物对环境因子的适应性进化[5–8],开展叶片解剖结构分析,对研究植物在极端生态环境下的适应策略具有重要意义。

同质园试验能最大程度地降低立地条件所造成的影响[9,10],进而说明不同种源植物叶片解剖结构特征表现出的差异是由于原种源地环境差异造成的植物遗传性基因变异。这与先前关于环境和遗传因素共同影响(包括气候梯度和种源的影响)叶片解剖结构特征的研究有所不同,这种同质园方法在局部适应性研究中具有重要意义[11],有助于我们更好地理解不同地理种源植物对全球变化的响应和适应机制。目前采用此法对不同种源内植物叶片解剖结构差异及其影响因素的相关分析较少。鉴于此,我们将来自华北和西北7个地区的蒙古莸种植在同一环境条件下,排除外界环境因素差异的作用,采用石蜡切片法对蒙古莸叶片解剖结构特征进行研究,以期进一步探讨植物的局域环境适应性进化,为今后蒙古莸的种质资源保护、遗传改良以及环境变异驱动的植物基因型差异对植物的作用机制提供理论依据。

2 数据集元数据简介

《7个种源蒙古莸叶片解剖结构及其地理环境影响因子数据集》[12]元数据信息,如名称、作者、地理区域、数据年代、数据集组成、数据共享政策等信息见表1

表1 《7个种源蒙古莸叶片解剖结构及其地理环境影响因子数据集》数据集元数据简表
条目描述
数据集名称7个种源蒙古莸叶片解剖结构及其地理环境影响因子数据集
数据集短名LeafAnatomicalStructure_CaryopterisMongholica
作者信息纪若璇 AAE-6059-2021, 北京林业大学, bljrx@bjfu.edu.cn
于笑 AAE-6050-2021, 北京林业大学, yuxiao@bjfu.edu.cn
常远, 北京林业大学, 1034530229@qq.com
沈超 0000-0001-5037-4922, 北京林业大学, shenchaonow@bjfu.edu.cn
白雪卡, 北京林业大学, baixueqia88@qq.com
夏新莉, 北京林业大学, xiaxl@bjfu.edu.cn
尹伟伦, 北京林业大学, yinwl@bjfu.edu.cn
刘超 AAE-6091-2021, 北京林业大学, liuchao1306@bjfu.edu.cn
地理区域38.88°N–45.73°N,105.72°E–116.79°E
数据年代2007–2017
数据格式.shp, .kmz, .xls, .jpg, .doc
数据量5.84 MB
数据集组成(1)种源地地理位置数据
(2)蒙古莸叶片横切面图片
(3)蒙古莸叶片解剖结构与地理环境影响因子分析数据,包括:
1) 种源地气候(年降水量、年平均气温、潜在蒸散量、海拔高度等)对蒙古莸叶片各个解剖结构参数的GLM分析数据
2) 蒙古莸叶片解剖性状相关性分析数据
3) 蒙古莸叶片解剖结构指标与种源地基本情况
基金项目生态环境部(2017ZX07101002);国家自然科学基金(32071734, 31770649, 31600484)
数据计算环境ArcGIS 10.2, Excel, R 3.2.2
出版与共享服务平台全球变化科学研究数据出版系统 http://www.geodoi.ac.cn
地址北京市朝阳区大屯路甲11号100101,中国科学院地理科学与资源研究所
数据共享政策全球变化科学研究数据出版系统的“数据”包括元数据(中英文)、通过《全球变化数据仓储电子杂志(中英文)》发表的实体数据和通过《全球变化数据学报(中英文)》发表的数据论文。其共享政策如下:(1)“数据”以最便利的方式通过互联网系统免费向全社会开放,用户免费浏览、免费下载;(2)最终用户使用“数据”需要按照引用格式在参考文献或适当的位置标注数据来源;(3)增值服务用户或以任何形式散发和传播(包括通过计算机服务器)“数据”的用户需要与《全球变化数据学报(中英文)》编辑部签署书面协议,获得许可;(4)摘取“数据”中的部分记录创作新数据的作者需要遵循10%引用原则,即从本数据集中摘取的数据记录少于新数据集总记录量的10%,同时需要对摘取的数据记录标注数据来源[13]
数据和论文检索系统DOI,DCI,CSCD,WDS/ISC,GEOSS,China GEOSS,Crossref

3 数据研发方法

3.1 研究地概况

同质园试验地点位于北京林业大学苗圃(40.01°N,116.34°E),属暖温带半湿润大陆性季风气候。年平均气温11.8 ℃,年最高气温37.5–42.6 ℃,最低气温–19.5––14.8 ℃,年降水量500–650 mm,多集中于夏季,占全年降水量的70%,而春旱情况严重[14]。试验地设置塑料遮雨棚以防止其他水分影响,雨天关闭,晴天打开[14]

3.2 试验材料

试验材料为蒙古莸1年生扦插苗。于2016年3月下旬,在内蒙古阿巴嘎旗、阿拉善左旗(阿左旗)、东乌珠穆沁旗(东乌旗)、鄂托克旗蒙西镇、凉城县、河北康保县、陕西神木市(种源地地理环境见表2)采取生长良好的蒙古莸植株,剪取长度约为20 cm左右的枝条茎段,放置于湿布中保存,立即带回北京。将不同种源地的扦插枝条经0.1%生根粉ABT1(中国林业科学研究院林业研究所)处理后扦插于450 mm×450 mm的花盆中。种植土为砂壤土(草炭土∶壤土∶沙子体积比为4∶4∶3),与蒙古莸的野外土壤成分相似,所有种源材料种植在同一土壤条件下。

表2 蒙古莸种源地基本情况
种源地纬度(°N)经度(°E)海拔(m)年降水量(mm)年平均气温(℃)生长季降水量(mm)生长季平均气温(℃)最冷月平均气温(℃)潜在蒸散量(mm)
内蒙古阿巴嘎旗43.90115.351,177224.032.9434.7117.74-19.85588.04
内蒙古阿左旗38.88105.721,670162.549.90315.5221.05-7.51712.81
内蒙古东乌旗45.73116.791,017200.722.46415.4617.76-20.52531.06
内蒙古蒙西40.08106.921,193118.008.92243.6420.95-9.58702.17
内蒙古凉城40.66112.301,429313.755.18611.4617.13-13.17580.80
陕西神木39.29110.331,209369.799.78729.8520.85-7.64706.91
河北康保41.99114.851,590279.723.82556.5316.78-15.22544.06

3.3 试验方法及算法原理

采用常规石蜡切片法[15]制作石蜡切片。将固定24 h以上的叶块,经70%、85%、95%、100%四个浓度酒精脱水。脱水时封盖,防止空气中水分渗入;后经二甲苯透明;将透明好的材料换入石蜡液浸透;浸透后的材料放入纸盒包埋(包埋所用石蜡熔点为56–57 ℃),待凝固后切片,切片厚度8–10 μm;番红-固绿双重染色,中性树胶封片,将其制作为永久切片[16]。在Leica DM2500显微镜(Leica, Wezlar, Germany)下观察,并用Leica LAS AF软件拍照记录。

用Image J软件中的直线工具测量叶片的上角质层厚度,上、下表皮细胞厚度,上、下栅栏组织厚度,叶厚度等叶解剖结构指标,并计算叶片结构紧密度、栅栏组织总厚度[17],每个指标重复测量5组数据。

叶片解剖结构指标:

叶片结构紧密度(Tight)= 栅栏组织厚度/叶片厚度 × 100%

栅栏组织总厚度(PT)= 上栅栏组织厚度(UPT)+ 下栅栏组织厚度(LPT

叶片解剖结构可塑性指标:

变异系数 = 标准差/算术平均值

可塑性指数 = (最大值-最小值)/最大值

根据蒙古莸不同种源地的经度与纬度,从ArcGIS 10.2的WorldCLIM全球高分辨率气候数据库[18]中选取各样地的月平均气候数据(2007–2016),并计算各气候指标大小。气候指标包括:年降水量(MAP)、生长季降水量(GSP)、降水年较差(ARP)、年平均气温(MAT)、生长季平均气温(GST)、气温年较差(ART)、潜在蒸散量(PE)。

3.4 数据分析

用方差分析多重比较法,比较不同种源蒙古莸叶片解剖结构参数差异。用Pearson相关性分析方法分析叶片解剖结构参数与种源地环境因子的相关关系。

结合一般线性模型(GLM) ANOVA分析种源地气候对蒙古莸叶解剖结构参数值的影响作用大小[19]。基于种源地环境因素与叶片解剖结构参数值的相关关系,选取MAPMATPE、海拔(Alt)、地理位置(Site)等5个种源地气候环境因子进入模型。

数据分析在Excel、R 3.2.2中完成。

3.5 技术路线

本数据集研发的主要步骤包括:

(1)基于前期资料收集与实地考察,采集7个种源蒙古莸植株,种植于北京林业大学苗圃;

(2)待植株生长一年半后取其叶片固定并切片,测量叶片的解剖结构指标,分析解剖性状间的相关性;

(3)分析原种源地地理环境信息与解剖结构参数间的相关性,应用GLM(一般线性模型)验证,形成7个种源蒙古莸叶片解剖结构及地理环境数据集。

技术路线如图2所示。

数据集研发技术路线图

图2 《7个种源蒙古莸叶片解剖结构及地理环境影响因子数据集》研发技术路线

下载

4 数据结果与验证

4.1 数据集组成

《7个种源蒙古莸叶片解剖结构及其地理环境影响因子数据集》由4部分文件组成。文件包括蒙古莸种源地的坐标数据、2007–2016年气候指标数据,叶片解剖结构示意图与参数信息,以及种源地气候与叶片解剖结构参数的GLM分析数据。数据文件存储为.shp、.kmz、.xls、.jpg和.doc格式。

4.2 数据结果

观察蒙古莸叶片解剖结构,7个种源的蒙古莸叶片均为等面叶,叶片厚度在192.34‒270.30 μm之间(表3);其内部结构可分为表皮、叶肉和叶脉3个基本部分(图3)。

表3 蒙古莸叶片解剖结构指标
解剖结构指标观测数目平均值标准偏差均值的95%置信区间最大值最小值变异系数可塑性指数
上角质层厚度 Cuticle(μm)706.871.666.05–7.7011.925.070.240.57
上表皮细胞厚度 UEC(μm)7018.535.3715.86–21.2026.469.870.290.63
下表皮细胞厚度 LEC(μm)7011.142.329.99–12.3017.648.740.210.50
上栅栏组织厚度 UPT(μm)70106.3517.9497.43–115.27146.6489.240.170.39
下栅栏组织厚度 LPT(μm)7076.1713.4469.48–82.85109.4259.570.180.46
栅栏组织总厚度 PT(μm)70182.5229.16168.02–197.02243.12156.050.160.36
叶厚度 LT(μm)70220.0225.71207.24–232.80270.30192.340.120.29
叶片结构紧密度 Tight(%)700.830.040.80–0.850.900.740.050.18
蒙古莸叶片横切面图

图3 蒙古莸叶片横切面(A, 叶肉结构。B, 主脉结构。C, 主脉及叶肉结构)
(Co,厚角组织;Cuticle,上角质层;LEC,下表皮细胞;LPT,下栅栏组织;MVP,主脉韧皮部;MVX,主脉木质部;St,孔下室;UEC,上表皮细胞;UPT,上栅栏组织)

下载

多数叶片解剖结构参数间存在明显的自相关,其中,UPTLPTPTLT以及Tight这5组参数两两之间都极显著正相关性(P < 0.01)(表4)。

表4 蒙古莸叶片解剖性状相关性分析
相关系数t上角质层厚度Cuticle(μm)上表皮细胞厚度UEC(μm)下表皮细胞厚度LEC(μm)上栅栏组织厚度UPT(μm)下栅栏组织厚度LPT(μm)栅栏组织总厚度PT(μm)叶厚度LT(μm)叶片结构紧密度Tight(%)
上角质层厚度 Cuticle(μm)0.2010.590**-0.246-0.434-0.351-0.260-0.432
上表皮细胞厚度 UEC(μm)0.2010.107-0.355-0.545*-0.470*-0.268-0.770*
下表皮细胞厚度 LEC(μm)0.590**0.107-0.203-0.436-0.326-0.259-0.391
上栅栏组织厚度 UPT(μm)-0.246-0.355-0.2030.721**0.948**0.945**0.683**
下栅栏组织厚度 LPT(μm)-0.434-0.545*-0.4360.721**0.904**0.825**0.806**
栅栏组织总厚度 PT(μm)-0.351-0.470*-0.3260.948**0.904**0.962**0.792**
叶厚度 LT(μm)-0.260-0.268-0.2590.945**0.825**0.962**0.595**
叶片结构紧密度 Tight(%)-0.432-0.770*-0.3910.683**0.806**0.792**0.595**

* P < 0.05; **P < 0.01.

4.3 数据结果验证

表5可得,蒙古莸叶片解剖结构指标与一个或多个种源地地理位置(经纬度)、温度、降水之间具有显著的相关性。

经度、纬度除与UEC呈正相关关系外,与UPTLPTPTLTTight这5个参数呈负相关关系(P < 0.05)。地理梯度上呈现自西向东、自南向北,UEC逐渐增大,UPTLPTPTLTTight逐渐减小的趋势。

表5 种源地环境因素与蒙古莸叶片解剖结构参数的相关系数
相关系数上角质层厚度Cuticle(μm)上表皮细胞厚度UEC(μm)下表皮细胞厚度LEC(μm)上栅栏组织厚度UPT(μm)下栅栏组织厚度LPT(μm)栅栏组织总厚度PT(μm)叶厚度LT(μm)叶片结构紧密度Tight(%)
纬度(°N)0.3280.638**0.451-0.526*-0.694**-0.643**-0.529*-0.724**
经度(°E)0.4370.704**0.366-0.787**-0.801**-0.854**-0.730**-0.886**
海拔(m)-0.124-0.034-0.3580.547*0.579*0.604**0.690**0.237
年降水量(mm)0.2610.074-0.103-0.671**-0.250-0.528*-0.555*-0.307
生长季降水量(mm)0.2440.085-0.110-0.689**-0.269-0.548*-0.575*-0.321
降水年较差(mm)0.3480.195-0.064-0.640**-0.236-0.503*-0.495*-0.365
年平均气温(℃)-0.446-0.803**-0.4210.534*0.735**0.668**0.495*0.860**
生长季平均气温 GST(℃)-0.476*-0.863*-0.2900.512*0.672**0.625**0.4240.880**
气温年较差 ART(℃)0.3380.503*0.529*-0.426-0.626**-0.551*-0.465-0.597**

* P < 0.05; ** P < 0.01; *** P < 0.001.

通过表6中ANOVA GLM分析可知,原种源地气候(温度、降水、蒸散)驱动的遗传变异对蒙古莸各解剖结构参数有显著作用,其解释度在34.09%–81.43%间,对LEC的解释度最小(34.09%),对Tight的解释度最大(81.43%)。

表6 种源地气候对蒙古莸叶片各个解剖结构参数的GLM分析
气候因子上角质层厚度Cuticle上表皮细胞厚度UEC
自由度df误差平方和均方误差F value解释度%SS显著性自由度df误差平方和均方误差F value解释度%SS显著性
年降水量10.0150.0152.5310.130.14713 × 10??3 × 10??0.070.120.804
年平均气温10.0310.0315.1820.770.049*10.2040.20434.8465.822.3 × 10??***
潜在蒸散量10.0060.0061.024.080.33910.0280.0284.879.200.055*
海拔10.0040.0040.682.740.43012.7832.7834.760.000.999
位置40.0390.0101.6326.200.24840.0240.0061.047.860.438
残差90.0540.006NA36.09NA90.0530.006NA17NA
总计170.150100170.310100
气候因子下表皮细胞厚度LEC叶厚度LT
自由度df误差平方和均方误差F value解释度%SS显著性自由度df误差平方和均方误差F value解释度%SS显著性
年降水量10.0010.0010.230.980.64610.0120.01285.5031.406.8 × 10??***
年平均气温10.0190.0193.8416.780.082*10.0090.00959.0521.693.0 × 10??***
潜在蒸散量10.0190.0193.7416.330.08510.0060.00639.0514.341.5 × 10??***
海拔19 × 10??9 × 10??0.020.080.89410.0070.00745.6216.768.3 × 10??***
位置40.0310.0081.5226.520.27640.0050.0018.5112.500.004**
残差90.0450.005NA39.30NA90.0011× 10??NA3.31NA
总计170.115100170.040100
气候因子上栅栏组织厚度UPT下栅栏组织厚度LPT
自由度df误差平方和均方误差F value解释度%SS显著性自由度df误差平方和均方误差F value解释度%SS显著性
年降水量10.0370.03793.7245.934.7 × 10??***10.0040.0047.474.780.023*
年平均气温10.0200.02050.1824.595.8 × 10??***10.0480.04890.2057.745.5 × 10??***
潜在蒸散量10.0050.00512.826.280.006**10.0050.0058.565.480.017*
海拔10.0070.00717.658.650.002**10.0070.00713.018.330.006**
位置40.0080.0025.1710.130.019*40.0150.0046.9917.910.008**
残差90.0044 × 10??NA4.41NA90.0055 × 10??NA5.76NA
总计170.081100170.084100
气候因子栅栏组织总厚度PT叶片结构紧密度Tight
自由度df误差平方和均方误差F value解释度%SS显著性自由度df误差平方和均方误差F value解释度%SS显著性
年降水量10.0190.019208.3927.381.6 × 10??***18 × 10??8 × 10??6.818.930.028*
年平均气温10.0310.031337.6644.371.9 × 10??***10.0070.00754.8071.894.1 × 10??***
潜在蒸散量10.0050.00551.636.785.2 × 10??***16 × 10??6 × 10??0.460.610.514
海拔10.0070.00773.499.661.3 × 10??***11 × 10??1 × 10??0.090.120.772
位置40.0070.00220.2010.621.6 × 10??***46 × 10??2 × 10??1.276.650.097*
残差98 × 10??9 × 10??NA1.18NA90.0011 × 10??NA11.81NA
总计170.070100170.009100

* P < 0.05; * P < 0.01; *** P < 0.001.

5 讨论与总结

本文研究了我国7个种源蒙古莸叶片解剖结构,分析了影响不同种源叶片解剖结构差异的重要因素,并发布了叶片解剖结构及其地理环境影响因子数据集。数据年代为2006–2017年,时间分辨率为天。利用同质园法最大程度降低了7个种源蒙古莸生长环境梯度的影响,结合多重比较、相关性分析、一般线性模型分析等方法,对7个种源蒙古莸叶片解剖结构及其影响因子进行定量比较,验证解剖结构参数与气候指标间的关系,最终得出结论,种源地气候差异驱动的遗传变异是引起不同种源叶片解剖结构差异的重要因素。本数据集丰富了局部适应性研究数据,帮助我们进一步了解不同种源植物应对全球变化的适应机制,可为全球变化、种质资源保护等研究领域提供基础数据。

作者分工:刘超、夏新莉、尹伟伦、纪若璇对数据集的开发做了总体设计;于笑、常远、沈超、白雪卡、纪若璇采集了植物材料;纪若璇处理了叶片相关数据并做了数据验证;刘超设计了模型和算法;纪若璇撰写了数据论文等。

利益冲突声明:本研究不存在研究者以及与公开研究成果有关的利益冲突。

参考文献